Metallicity Gradients in the Intracluster Gas of Abell 496 

Renato A. Dupke^'^ & Raymond E. White III^'^ 
^Department of Physics & Astronomy, University of Alabama, Tuscaloosa, AL 35487-0324 

^Department of Astronomy, University of Michigan, Ann Arbor, MI 48109-1090 
^Laboratory for High Energy Astrophysics, Code 662, NASA/GSFC, Greenbelt, MD 20771 
O" 

o 

O . ABSTRACT 
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r.^ I Analysis of spatially resolved ASCA spectra of the intracluster gas in Abell 

^ ' 496 confirms there are mild metal abundance enhancements near the center, 
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CSJ I as previously found by White et al. (1994) in a joint analysis of Ginga LAC 
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and Einstein SSS spectra. Simultaneous analysis of spectra from all ASCA 



*> \ instruments (SIS + CIS) shows that the iron abundance is 0.36 ± 0.03 solar 

3 — 12' from the center of the cluster and rises ~50% to 0.53 ± 0.04 solar within 
the central 2'. The F-test shows that this abundance gradient is significant at 

CN ■ 6 6 

O I the >99.99% level. Nickel and sulfur abundances are also centrally enhanced. 



We use a variety of elemental abundance ratios to assess the relative contribution 
of SN la and SN II to the metal enrichment of the intracluster gas. We find 



Qh| spatial gradients in several abundance ratios, indicating that the fraction of 

O ' iron from SN la increases toward the cluster center, with SN la accounting 

5-1 ' 

^ ' for ~50% of the iron mass 3 — 12' from the center and ~70% within 2'. The 

f^ ■ increased proportion of SN la ejecta at the center is such that the central 

iron abundance enhancement can be attributed wholly to SN la; we find no 
/\ • significant gradient in SN II ejecta. These spatial gradients in the proportion 

c^ ■ of SN la/II ejecta imply that the dominant metal enrichment mechanism near 

the center is different than in the outer parts of the cluster. We show that the 
central abundance enhancement is unlikely to be due to ram pressure stripping 
of gas from cluster galaxies, or to secularly accumulated stellar mass loss within 
the central cD. We suggest that the additional SN la ejecta near the center is 
the vestige of a secondary SN la-driven wind from the cD (following a more 
energetic protogalactic SN Il-driven wind phase), which was partially smothered 
in the cD due to its location at the cluster center. 



Subject headings: galaxies: abundances — galaxies: clusters: individual (Abell 
496) — intergalactic medium — X-rays: galaxies 
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1. Introduction 

While tlie metals observed in intracluster gas clearly originate from stars, it remains 
controversial how the metals got from stars into the intracluster gas. The two global metal 
enrichment mechanisms considered to be most likely are protogalactic winds from early-type 
galaxies (Larson & Dinerstein 1975) and ram pressure stripping of gas from galaxies (Gunn 
& Gott 1972). Early Einstein FPCS spectroscopy (Canizares et al. 1982) and more recent 
ASCA spectroscopy (Mushotzky & Loewenstein 1997; Mushotzky et al. 1996) showed that 
global intracluster metal abundances are consistent with ejecta from Type II supernovae, 
which supports the protogalactic wind model. White (1991) showed that the specific 
energy of intracluster gas is greater than that of cluster galaxies, which also suggests that 
protogalactic winds injected significant amounts of energy and metals into intracluster gas. 
However, theoretical uncertainties about the elemental yields from Type II supernovae make 
it difficult to determine confidently the relative proportion of iron from SN II and SN la in 
intracluster gas (Gibson, Loewenstein & Mushotzky 1997). This uncertainty has allowed 
others to conclude that as much as ~50% of the iron in intracluster gas comes from SN la 
(Ishimaru & Arimoto 1997; Nagataki & Sato 1998). The possible presence of such large 
quantities of iron from SN la is problematic. Is ram pressure stripping so effective that it 
contaminates the outer parts of clusters nearly as effectively as the central regions? Or is it 
ejecta from a secondary SN la-driven wind phase in ellipticals? Clues about the dominant 
enrichment mechanism(s) may be found in the detailed spatial distribution of elements in 
intracluster gas. 

An increasing number of galaxy clusters are being found with centrally enhanced metal 
abundances in their intracluster gas. Ginga observations of the Virgo cluster showed that its 
iron abundance declines from ~0.5 solar at the center to ~0. 1-0.2 solar 3° away (Koyama, 
Takano & Tawara 1991). White et al. (1994) found central abundance enhancements in 
Abell 496 and Abell 2142 in joint analyses of Ginga LAC and Einstein SSS spectra. ASCA 
observations of the Centaurus cluster show that its iron abundance declines from ~ solar 
at the center to ~0.3 solar 15' away (Fukazawa et al. 1994). The Perseus cluster also has 
an abundance gradient near its center (Ulmer et al. 1987; Ponman et al. 1990; Kowalski 
et al. 1993; Arnaud, et al. 1994). More recently, central abundance enhancements were 
found in ASCA data for Hydra A (Ikebe et al. 1997), AWM 7 (Ezawa et al. 1997; Xu et al. 
1997), Abell 2199 and Abell 3571 (Dupke 1998, Dupke & White 2000), as well as in ROSAT 
(Pislar et al. 1997) and BeppoSax (Irwin & Bregman 1999) data for Abell 85. The presence 
of these central abundance enhancements is poorly correlated with global cluster properties. 

In this paper we analyze spatially resolved ASCA spectra of Abell 496 and confirm that 
it has centrally enhanced metal abundances. Abell 496 is a Bautz-Morgan Type I cluster 
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with an optical redshift oi z = 0.0328. Adopting a Hubble constant of 50 km s~^ Mpc^^ and 
go = 0.5, its luminosity distance is 197 h'^Q Mpc and 1' = 57 h^Q kpc. The central cD (MCG 
-02-12-039) has a total B magnitude of Bt = 13.42 (Valentijn 1983) and an optical effective 
radius of rcfr > 49 h'^Q kpc (Schombert 1986). Neither the projected galaxy distribution nor 
the galaxy velocity distribution in the cluster shows signs of significant substructure, so the 
cluster appears to be dynamically relaxed (Bird 1993; Zabludoff, Huchra & Geller 1990). 
Heckman (1981) was the first to suggest that Abell 496 contained a cooling flow, after he 
detected cool Ha-emitting gas in its central cD galaxy; subsequent optical observations 
found extended Ha emission (Cowie et al. 1983). Nulsen et al. (1982) found a soft X-ray 
component in Einstein SSS spectra of this cluster and estimated a cooling accretion rate 
of ~200 Mq jt~^, which is consistent with later analyses (Mushotzky 1984; Mushotzky 
& Szymkowiak 1988; Canizares, Markert & Donahue 1988; Thomas, Fabian & Nulsen 
1987; White et al. 1994). In the course of a joint analysis of Einstein SSS and Ginga 
LAC spectra. White et al. (1994) found a central abundance enhancement in Abell 496; 
the differing fields of view of these two instruments allowed a coarsely spatially-resolved 
analysis. The ASCA observations of Abell 496 have been previously analyzed by Mushotzky 
(1995) and Mushotzky et al. (1996), who found no evidence of an abundance gradient 
within the central 1 Mpc; however, the first paper analyzed data from only one of the four 
ASCA spectrometers, while the second considered data from all four ASCA spectrometers, 
but only beyond 3' from the center (in order to avoid the spectral influence of the central 
cooling flow). We include the central cooling flow region in our analysis of data from all 
four ASCA spectrometers. 



2. Data Reduction &; Analysis 

Abell 496 was observed for 40 ksec by ASCA on 20-21 September 1993. ASCA 
carries four large-area X-ray telescopes, each with its own detector: two Gas Imaging 
Spectrometers (GIS) and two Solid-State Imaging Spectrometers (SIS). Each GIS has a 50' 
diameter circular field of view and a usable energy range of 0.7-12 keV; each SIS has a 
22' square field of view and a usable energy range of 0.4-10 keV. We selected data taken 
with high and medium bit rates, with cosmic ray rigidity values > 6 GeV/c, with elevation 
angles from the bright Earth of > 20° and from the Earth's limb of > 5° (GIS) or 10° (SIS); 
we also excluded times when the satellite was affected by the South Atlantic Anomaly. Rise 
time rejection of particle events was performed on GIS data, and hot and flickering pixels 
were removed from SIS data. The resulting effective exposure times for each instrument are 
shown in Table 1. Since the cluster flUs the spectrometers' fields of view, we estimated the 
background from blank sky files provided by the ASCA Guest Observer Facility. 
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We used XSPEC v9.0 (Arnaud 1996) software to analyze the SIS and GIS spectra 
separately and jointly. We fit spectra using the mekal and vmekal thermal emission models, 
which are based on the emissivity calculations of Mewe & Kaastra (cf. Mewe, Gronenschild 
& van den Oord 1985; Mewe, Lemen & van den Oord 1986; Kaastra 1992), with Fe L 
calculations by Liedahl, Osterheld & Goldstein (1995). Abundances are measured relative 
to the solar photospheric values of Anders & Grevesse (1989), in which Fe/H=4.68 x 10~^ 
by number. Galactic photoelectric absorption was incorporated using the wabs model 
(Morrison & McCammon 1983); the Galactic column of absorbing material in this line of 
sight is Nh = 4.58 x 10^° cm^^ (Dickey & Lockman 1990; HEASARCllR software). Spectral 
channels were grouped to have at least 25 counts/channel. Energy ranges were restricted to 
0.8-10 keV for the GIS and 0.4-10 keV for the SIS. The minimum projected sizes for regions 
of spectral extraction (circular radii or annular widths) were typically 3' for the GIS and 2' 
for the SIS. To maximize the statistical significance of any gradients, we also assessed larger 
regions in the outer parts. A central region of GIS data with a projected radius of 2' was 
also analyzed for closer comparison with the SIS. The intracluster gas temperature in Abell 
496 is cool enough (~4 keV) that the energy dependence of ASCA's point spread function 
does not affect our results significantly. Since our results for spectral fits to the individual 
instruments are consistent with our results from the joint analysis of all instruments, we 
will discuss only the joint analysis. 



3. Results 

3.1. Temperature and Abundance Profiles 

We jointly fitted thermal models to spectra from all four ASCA instruments. The 
fitting normalizations for the data from each instrument were allowed to vary independently, 
in order to compensate for small calibration and spatial extraction differences between 
the detectors; the normalizations differ by < 10% in practice. Figure la shows the GIS 
data with the best fits in both the inner and outer projected regions of the cluster, while 
Figure lb shows the same fits to the SIS data from these regions. The resulting fits had 
X^ ~ 1 and the temperature and abundance distributions are shown in Figures 2a & b 
and Table 2. The temperature rises from 3.24to;o6 keV within 2' to 4.40lo;i3 keV beyond 
5'. The central abundance is 0.53to;o4 solar, falling to 0.36to;Q3 in the outer 3-12' region. 
In this combined data set, the central abundance enhancement within 3' is significant at a 
confidence level of ~99%. Confidence contours for the abundances in the two innermost 
regions are shown in Figure 3. We also used the F-test to assess the significance of the 
central abundance enhancement. We simultaneously fitted the spectra from the inner 



-5- 



and outer regions, allowing their normalizations, temperatures and abundances to vary 
independently. We then refit the spectra with the abundances from the two regions tied 
together. The x^ of these latter fits were larger by > 26 for an increase of only one degree 
of freedom. The F-test indicates that the central abundance enhancement is significant at 
a level of >99.99%. 

Since there is a moderate cooling flow at the center of Abell 496, we tested whether the 
abundance gradient described above is an artifact of our choice of spectral model. We added 
a cooling flow component to the mekal thermal emission model in the central region. The 
cooling flow spectral model cf low in XSPEC is characterized by maximum and minimum 
temperatures, an abundance, a slope which parameterizes the temperature distribution 
of emission measures, and a normalization, which is simply the cooling accretion rate. 
We adopted the emission measure temperature distribution that corresponds to isobaric 
cooling flows (zero slope). We tied the maximum temperature of the cooling flow to the 
temperature of the thermal component, and we fixed the minimum temperature at 0.1 
keV. We applied a single (but variable) global absorption to both spectral components and 
associated an additional, intrinsic absorption component with the cooling flow, placing 
it at the redshift of the cluster. The addition of the cooling flow component does not 
significantly affect our results for the central region: the central abundance enhancement 
remains significant at a confidence level of >90%. The two component fits at the center are 
slightly worse than the isothermal model fits above, but they still have x^ ~ 1 (see Table 
2). In order to apply the F-test to these cooling fiow model fits, we also simultaneously fit 
the spectra from inner and outer regions, first allowing the abundances in the respective 
regions to vary independently, then tying the abundances together. The F-test implies that 
the abundance gradient is significant at the >99.99% level. Since the spectra we analyze are 
from cylindrical projections of emission measure through various lines of sight through the 
cluster, the true spatial abundance enhancement at the center will be somewhat stronger 
than we observe. 



3.2. Individual Elemental Abundances & Abundance Ratios 

We also determined the abundances of individual elements using the vmekal spectral 
model in XSPEC. A similar analysis for the outer regions (> 3') of Abell 496 was done by 
Mushotzky et al. (1996), whose individual elemental abundance measures are consistent 
with our results at the 90% confidence level. In our spectral model fits, the He abundance 
was fixed at the solar value, while C and N were fixed at 0.3 solar (since ASCA is rather 
insensitive to C and N and the derived abundances of other elements are not affected by the 
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particular choice for C and N abundances). Our observed abundances are shown in Table 3 
for various projected spatial regions. 

Table 3 shows that the iron abundance is best determined and increases ~50% from 
0.36 ± 0.03 solar in the outer parts (> 3') to a central value of 0.53 ± 0.05 solar. The 
sulfur abundance also shows a significant gradient, rising from 0.20±0.16 solar beyond 3' 
to 0.58±0.20 solar at the center. The silicon abundance is ~0.8±0.2 solar, showing no 
significant gradient within its 90% confidence limits. The best fitting neon abundance is 
also nearly solar, showing no significant gradient, while the best fitting nickel abundance is 
~2.5 times solar at the center, with a marginally significant decline to solar in the outer 
parts. Oxygen, magnesium, argon and calcium are poorly constrained. 

Theoretical numerical models of supernova yields predict the following elemental ratios 
relative to solar values: for SN la, the W7 model of Nomoto, Thielemann & Yokoi (1984), 
as updated in Nomoto et al. (1997a), gives 

O ^ Mg ^ 0.035 Fe, 

Ne ^ 0.006 Fe, 
Si ^ S ^ Ar ^ Ca ^ 0.5 Fe, 
Ni ^ 4.8 Fe; 
while for SN II, Nomoto et al. (1997b) find 

O ^ Mg ^ Si ^ 3.7 Fe, 

Ne ^ S ^ 2.5 Fe, 

Ar ?« Ca ?« Ni ^ 1.7 Fe, 

after integrating their yields over a Salpeter mass function with upper and lower mass limits 
of 10 and 50 Mq, respectively. 

Various observed abundance ratios within inner and outer projected spatial regions 
are shown in Table 4, along with the theoretical expectations for SN la and SN II ejecta; 
the errors associated with the observed abundance ratios are the propagated la errors. 
Note that several abundance ratios lie significantly outside their theoretical ranges: Ne/Si 
(0 — 3'), Ne/S (3 — 12') and Si/S (3 — 12'). Three abundance ratios show marginally 
significant gradients: Si/S, Si/Ni, and S/Fe. 
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4. Distinguishing the Relative Contributions from SN la & SN II 

We use the observed abundance ratios (normalized by their solar values) in Table 4 
to estimate the relative contribution of SN la and SN II to the metal enrichment of the 
intracluster gas. Such estimates are complicated by uncertainties in both the observations 
and the theoretical yields. The yield relations above show that the most discriminatory 
abundance ratios are the ones involving oxygen, magnesium and neon, since their ratios to 
iron are 2-3 orders of magnitude smaller for SN la than for SN II; of these, magnesium 
is poorly determined, which leaves oxygen and neon. Despite their large fractional 
uncertainties, the fact that the observed 0/Fe and Ne/Fe ratios are of order unity (see 
Table 4) clearly indicates the presence of SN II ejecta — these ratios are predicted to be 
less than a few percent for pure SN la ejecta. The 90% errors in these abundance ratios are 
not so large that they can be consistent with SN la ejecta alone. On the other hand, the 
observed 0/Fe and Ne/Fe ratios are about half the values predicted for pure SN II ejecta, 
which may indicate dilution by SN la ejecta. The SN la iron mass fractions indicated by 
these abundance ratios are indicated in Table 4. However, the theoretical iron yields from 
SN II are uncertain by a factor of ~2 (Woosley & Weaver 1995; Arnett 1996; Gibson et al. 
1997), so this systematic uncertainty obscures the relative contribution from SN la and SN 
II. 

The abundance ratios with the smallest fractional errors are Si/Fe, S/Fe and Ni/Fe. 
The Si/Fe and Ni/Fe ratios in Table 4 indicate a roughly comparable mix of SN II and 
SN la ejecta in the outer parts, while the SN la/II mix indicated by the S/Fe ratio is 
inconsistent with those derived from the Si/Fe and Ni/Fe ratios (we show below that 
elemental ratios involving sulfur are problematic). The associated SN la iron mass fractions 
are listed in Table 4 for these ratios as well. These three abundance ratios also involve iron, 
the production of which in SN II models is uncertain by a factor of ~2, so we next consider 
SN la/II discriminators that do not involve iron. 

Two well-determined ratios independent of iron are Si/Ni and Si/S. The SN la/II ratio 
derived from Si/Ni is consistent with the values derived above (see Table 4), indicating 
that SN la contribute 74% (66%) of the iron mass in the inner — 2' (0 — 3') region and 
48% in the outer 3 — 12' region. However, the best-fit value of the Si/S ratio is outside 
the theoretical boundaries in the outer parts of Abell 496, although its errors are large 
enough to be consistent with the expectation for SN II ejecta. Inspection of Table 4 shows 
that most of the best fit values of the other ratios involving sulfur in the outer parts are 
also systematically outside the theoretical range. We conclude, despite the large fractional 
errors, that sulfur is likely to be overproduced in the SN II models we have adopted. To be 
consistent with the results from the other elemental ratios above, sulfur production in SNe 



II should be reduced by a factor of ~2-4 relative to the models of Nomoto et al. (1997b). 
This overproduction of sulfur in other SN II models was noted previously by Mushotzky 
et al. (1996), Loewenstein & Mushotzky (1996) and Gibson et al. (1997). Recently, 
Nagataki & Sato (1998) found that this sulfur discrepancy was reduced when they used the 
theoretical SN II yields of Nagataki et al. (1998), who explored the affects of asymmetric 
SN II explosions. 

Figure 4 summarizes our estimates of the iron mass fraction from SN la, as derived 
from the variety of elemental abundance ratios described above; values for the inner — 2' 
region are indicated by filled circles, while empty circles correspond to 3 — 12'. The filled 
and empty symbols are obviously segregated, despite the large individual errors, indicating 
that the outer region is more dominated by SN II ejecta and the central region is dominated 
by SN la ejecta. The consistency of the results for ratios not involving sulfur is particularly 
noteworthy; in Figure 4 the average theoretical yield of sulfur from SN II has been reduced 
by a factor of 3.8 to bring the estimates from sulfur ratios more in line with the other 
estimates. 

For our consensus estimate of the iron mass fraction due to SN la, we average the SN 
la iron mass fractions for five of the seven abundance ratios discussed above: 0/Fe, Ne/Fe, 
Si/Fe, Ni/Fe and Si/Ni (the two ratios involving sulfur are excluded). Our best estimates 
for the SN la iron mass fraction are 70 ± 5% within — 2', 64 ± 5% within — 3', and 50 ± 6% 
within 3 — 12' (these values are denoted as /x in Figure 4). Although the exact proportions of 
SN la and SN II ejecta are sensitive to our adopted theoretical SN II yields, this constitutes 
evidence for an increase in the proportion of iron from SN la at the center relative to the 
outer parts. For our adopted yields, the proportion of iron from SN la is ~50% larger at 
the center than the outer parts. This increased proportion of SN la ejecta is such that the 
central iron abundance enhancement can be attributed wholly to SN la. Combining our 
results for the gradients in both the iron abundance and the relative proportion of SN la 
and II, we find that the iron abundance due to SN II is Apcn = 0.18 ± 0.06 solar in the 
outer parts (3 — 12') and 0.16 ± 0.06 solar near the center (0 — 2'), exhibiting no significant 
gradient. The iron abundance due to SN la is A^e^^ = 0.18 ± 0.06 solar in the outer parts, 
increasing a factor of ~2 to 0.38 ± 0.07 solar near the center (see Table 5). 

While our estimate of the relative proportion of SN la/II ejecta depends upon our 
adopted SN yield models, our qualitative conclusion that there is a large fraction of iron 
from SN la is robust. Fukazawa et al. (1998) found in ASCA spectra of 40 clusters that 
Si/Fe ratios are lower in cool clusters than in hot clusters, indicating that the proportion of 
SN la ejecta is larger in cool clusters than in hot clusters. Davis, Mulchaey & Mushotzky 
(1999) recently showed that this trend continues in galaxy groups dominated by early-type 
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galaxies. Given that both SN la and SN II produce significant amounts of Si and Fe 
(although in different proportions), there must be roughly comparable amounts of ejecta 
from both SN Types in typical clusters to produce the Si/Fe - kT trend. Fukazawa et al. 
(1998) conclude that the SN la/II mixture in clusters cannot be as ambiguous as suggested 
by Gibson et al. (1997). 



5. Distinguishing between Possible Metal Enrichment Mechanisms 

We showed above that iron, nickel and sulfur abundances are centrally enhanced 
in Abell 496 and that gradients in various elemental ratios indicated that these central 
abundance enhancements are largely due to SN la ejecta. Our more model dependent result 
is that ~50% of the iron 3 — 12' from the cluster center comes from SN la. The existence of 
spatial gradients in abundance ratios implies that the dominant metal injection mechanism 
near the cluster center must be different than in the outer parts. We will first distinguish 
between several possible mechanisms for producing central abundance enhancements in the 
cluster. Then we will assess the relative roles of winds and ram pressure stripping as metal 
enrichment mechanisms for the bulk of the cluster. 



5.1. Mechanisms for Creating Central Abundance Enhancements 

There are several mechanisms that may cause central abundance enhancements in 
intracluster gas: 1) ram pressure stripping of the metal-rich gas in cluster galaxies by 
intracluster gas is more effective at the center, where the intracluster gas density is highest 
(Nepveu 1981); 2) secular mass loss from the stars in central dominant galaxies may 
accumulate near the cluster center (White et al. 1994); 3) if galaxies blew winds which 
were not thoroughly mixed in the intracluster gas, metal abundances would decline outward 
from the center, since the luminosity density of metal- injecting galaxies falls more rapidly 
with radius than the intracluster gas density (Koyama et al. 1991; White et al. 1994); 
4) even if cluster galaxies blew winds which were generally well-mixed in the intracluster 
gas, a central dominant galaxy's wind may be at least partially suppressed, by virtue of 
its location at the bottom of the cluster's gravitational potential and in the midst of the 
highest intracluster gas density (White et al. 1994). Although each of these mechanisms 
can produce abundance gradients, no individual mechanism can produce spatial gradients 
in abundance ratios. The fact that we observe gradients in abundance ratios implies that 
the dominant metal enrichment mechanism changes spatially. 
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There are several reasons why ram pressure stripping is unlikely to be the source of 
SN la iron at the center of the cluster. First, the gaseous abundances measured in most 
early-type galaxies by ASCA (Loewenstein et al. 1994; Matsumoto et al. 1997) and ROSAT 
(Davis & White 1996) are 0.2-0.4 solar, significantly less than the 0.5-0.6 solar abundance 
observed at the cluster center. Only the most luminous ellipticals, which also tend to be 
at the centers of galaxy clusters or groups, are observed to have gaseous abundances of 
0.5-1 solar. Second, if ram pressure stripping is the primary source of SN la material at the 
center of the cluster, the cD is the one galaxy in the cluster which should not be stripped. 
Therefore, the cD should exhibit its accumulated history of SN la ejecta, in addition to any 
ejecta stripped from other galaxies as they passed near the cluster center. The cD should 
then have more SN la iron than expected for its luminosity. However, we show in Appendix 
A that the SN la iron mass to luminosity ratio {M-pci^/Lbj^,^^,) in the vicinity of the cD is 
no greater than in the rest of the cluster (see Table 5). This indicates that the bulk of the 
SN la ejecta from the cD has been retained in its vicinity, but has not been supplemented 
by ejecta from other galaxies. Third, gaseous abundances in ellipticals tend to decline 
outward from their centers, as observed in NGC 4636 (Matsushita et al. 1997) and other 
early-type galaxies (Matsushita 1997), since the mass-losing stars in ellipticals exhibit 
such metallicity gradients; thus, most stripped gas will have even lower abundances than 
indicated by the global X-ray measures of ellipticals (since global measures are weighted by 
centrally concentrated emission measures in ellipticals), which are already too low. Finally, 
since the efficiency of ram pressure stripping depends on the intracluster gas density, which 
declines strongly with radius, the abundances of SN la ejecta should also decline strongly 
with radius, which is not observed; we see only a factor of two decline in the iron abundance 
from SN la (see Table 5). We conclude that ram pressure stripping is not the cause of the 
central abundance enhancements in Abell 496. If ram pressure stripping is effective in the 
cluster, it would act to dilute the central abundance enhancement. 

To assess whether secularly accumulated stellar mass loss in the cD could cause the 
observed central abundance enhancement in Abell 496, we compare the cD to a giant 
elliptical, which is not in a rich cluster: NGC 4636 is one of the most X-ray luminous 
ellipticals and may be at the center of its own small group (Matsushita et al. 1998). Thus, 
NGC 4636 is not in an environment where it is likely to have been stripped. If giant 
(unstripped) ellipticals secularly accumulate their stellar mass loss (after a putative wind 
phase), NGC 4636 should have the same SN la iron mass to light ratio as the cD in Abell 
496. Instead, we find that the SN la Fe mass to light ratio in the central region of Abell 
496 is at least ~ 20 times higher than in NGC 4636 (see Appendix A). We conclude that 
the central abundance enhancements in Abell 496 are not likely to be due to secularly 
accumulated stellar mass loss in the cD. 
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If early-type galaxies blew winds which were only locally mixed in the intracluster gas, 
we would expect to observe an abundance gradient A{r) which is proportional to the ratio 
of the luminosity density of wind-blowing galaxies (^e/so) to the intracluster gas density, 
i.e. A{r) oc iE/so{r) / Pgas{r) . Equivalently, A{r) oc LE/so(^)/Mgas(r), where LE/so(r) and 
^gas(^) are the luminosity of early-type galaxies and the intracluster gas mass enclosed 
in spherical shells about the center. We derived the luminosity distribution of early-type 
galaxies in Abell 496 from the galaxy morphology data of Dressier (1980) and the gas 
density distribution from Einstein IPC data (see Appendix A). In Table 5 we show that 
the ratio -^E/so/^gas declines a factor of ~ 4 from — 2' to 3 — 12', while the observed 
iron abundance declines only 30% (and the SN la iron abundance drops by half). Thus, 
the observed abundance gradient is too shallow to be caused by poorly mixed winds from 
early-type galaxies in the cluster. 

We propose instead that the central gradients in abundances and abundance ratios in 
Abell 496 result from a partially suppressed SN la-driven wind from the cD. This would be 
a secondary wind phase, following a more vigorous, unsuppressed protogalactic wind driven 
by SN II. Such a secondary wind phase in noncentral early-type galaxies can generate 
the SN la enrichment seen in the bulk of the intracluster gas, as we suggest in the next 
subsection. As we showed above, the SN II iron mass to light ratio in the vicinity of the cD 
and the lack of a central enhancement in SN II ejecta indicate that it lost the bulk of its 
SN II ejecta. However, a weaker SN la-driven wind would be more readily suppressed at 
center of the cluster, due to the depth of the gravitational potential and the high ambient 
intracluster gas density. SN la-driven winds would be less vigorous than the initial SN 
Il-driven winds, since SN la inject ~10 times less energy per unit iron mass than SN II; the 
observations indicate that comparable amounts of iron came from SN la and SN II, so SN 
II have injected ~10 times more energy into the intracluster gas than SN la. 

If suppressed winds from central dominant galaxies are the cause of central abundance 
enhancements in other clusters, the prevalence of such enhancements in cooler clusters 
may be related to their cooling flow properties. Cool clusters tend to have cooling flows 
with smaller accretion rates than hot clusters, so the history of prior metal ejection from 
the central cD may be more likely to survive in cool clusters. For a given cD optical 
luminosity, the metal ejection is more likely to extend beyond the cooling flow region (since 
the central intracluster gas density is smaller in cool clusters than hot clusters), and the 
inward advection of the cooling flows would be less destructive to preexisting abundance 
gradients than in hotter clusters with higher accretion rates. 
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5.2. Global Enrichment Mechanisms 

As mentioned in the introduction, the two metal enrichment mechanisms usually 
considered for the bulk of intracluster gas are protogalactic winds from early-type galaxies 
and ram pressure stripping of gas from galaxies in the cluster. The yields of the SN II 
models that we adopt (Nomoto et al. 1997b) lead us to conclude that nearly 50% of the 
iron in Abell 496 comes from SN la. Similar conclusions were reached for other clusters by 
Ishimaru & Arimoto (1997) and Nagataki & Sato (1998), who used different theoretical 
models for SN II. Slightly more than half of the cluster iron comes from SN II, which can 
be readily attributed to protogalactic wind enrichment. However, the quantity and spatial 
extent of the iron from SN la is problematic: is ram pressure stripping so effective that it 
contaminates the outer parts of clusters nearly as effectively as the central regions? Or was 
there a (secondary) galactic wind phase driven by SN la? 

If ram pressure stripping is effective in the outer parts of the cluster, it should be 
even more effective at the center, where the intracluster gas density is highest. However, 
we showed in the previous subsection that ram pressure stripping cannot account for the 
central concentration of metals. Renzini et al. (1993) have also argued strongly against ram 
pressure stripping being very significant in clusters, citing the lack of a strong metallicity 
trend with cluster temperature: hot clusters have higher velocity dispersions and tend to 
have much higher gas densities than cool clusters, so ram pressure stripping should be much 
more effective in hot clusters than cool clusters. The lack of a strong metallicity trend with 
cluster temperature implies that ram pressure is not the major source of intracluster metals. 

We propose instead that the bulk of intracluster gas is contaminated by two phases of 
winds from early-type galaxies: an initial SN Il-driven protogalactic wind phase, followed 
by a secondary, less vigorous SN la-driven wind phase. As mentioned in the previous 
subsection, secondary SN la-driven winds would be ~ 10 times less energetic than the 
initial SN Il-driven protogalactic winds. Fukazawa et al. (1998) invoked SN Il-driven 
protogalactic winds to account for their discovery that the proportion of SN la ejecta is 
higher in cool clusters than in hot clusters. They suggested that protogalactic winds were 
energetic enough that SN Il-enriched material was able to escape cool clusters, which have 
shallower gravitational potentials than hot clusters. Less vigorous secondary SN la-driven 
winds would allow SN la-enriched material to escape most galaxies, but not clusters. 

The two phase wind scenario we are advocating has not been explicitly modeled to 
date. Some previous evolutionary models for the gas in ellipticals have presumed the 
existence of an initial SN Il-driven wind, without modeling it, and concentrated on a later 
SN la-drive wind phase (Loewenstein & Mathews 1991; Ciotti et al. 1991). Some of these 
models assume that the bulk of intracluster iron comes from SN la. Others investigators 
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have assumed that the bulk of intracluster iron comes from SN II and model only SN 
Il-driven protogalactic winds in detail (Larson & Dinerstein 1975; David, Forman, & Jones 
1991). However, in all recent investigations, present epoch SN la rates were adopted 
which lead to huge overpredictions of the current gaseous abundances in ellipticals: iron 
abundances were theoretically predicted to be 3-5 times solar, while ASCA and ROSAT 
observations find abundances to be 0.1-1 solar. 

In parlance similar to that of Renzini et al. (1993), we are proposing a wind-outflow- 
inflow (WOI) model in which the "wind" is driven by SN II, the less vigorous "outflow" 
is driven by SN la, and the subsequent inflow experiences much less contamination by SN 
la than in previous modeling. Generating SN la-enriched outflows that can contaminate 
intracluster gas, but leave current abundances subsolar in elliptical atmospheres, requires 
rather different SN la rate evolution than in previous models. Prior models of SN la-driven 
winds tended to inject roughly the right amount of iron in clusters (to within a factor of 
~2 or so). However, the SN la rate must decline much faster than in previous models if the 
current SN la rate is to be as low as ^ 0.03 SNU (Loewenstein & Mushotzky 1997), in order 
to match the low iron abundances in elliptical atmospheres. This rate is 3-10 times smaller 
than previously adopted and is ~ 4 times smaller than the most recent optical estimate 
of the current SN la rate in ellipticals (0.13 SNU; Capellaro et al. 1997). To generate the 
amount of SN la observed in intracluster gas, the SN la rate at earlier times must be much 
larger than previously modeled, to compensate for the much lower current rate. These 
heuristic constraints on the evolution of the SN la rate are not yet theoretically motivated. 



6. Summary 

We have carried out a detailed analysis of the distribution of elemental abundances in 
the intracluster gas of Abell 496. Our main results which are independent of our choice of 
supernovae yield models are: 

1. The hot gas of Abell 496 has significant abundance gradients: the iron abundance is 
0.361o;q| solar 3-12' from the center, rising ~50% to 0.53lo;o4 solar within 2'; nickel 
and sulfur also have significant central concentrations. 

2. There are spatial gradients in elemental abundance ratios in this cluster; a variety 
of abundance ratios individually and collectively indicate that SN la ejecta is more 
dominant in the center than in the outer parts. 

3. We find no significant gradient in SN II ejecta. 
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4. Ram pressure stripping is unlikely to generate the observed central abundance 
enhancements for several reasons, including the fact that gaseous abundances observed 
in elliptical atmospheres tend to be substantially less than the abundances observed 
in the intracluster gas near the center. 

5. Two stage galactic winds, consisting of SN Il-driven protogalactic winds followed by 
less energetic SN la-driven outflows, are proposed to generate comparable levels (by 
mass) of iron contamination from SN la and SN II in the intracluster gas of Abell 496. 

6. Since the secondary SN la-driven wind is ~ 10 times less energetic than the SN 
Il-driven wind, it is more likely to be smothered due to the cD being at the bottom 
of the cluster's gravitational potential and in the midst of the highest intracluster 
gas density; Such a smothered wind may generate the observed central SN la iron 
abundance enhancement. 

Results which are more dependent upon our particular choice of supernovae yield 
models include: 

1. SN la account for ~50% of the iron mass 3 — 12' from the center of the cluster and 
~70% of the iron mass within 2'. 

2. The central iron abundance enhancement can be attributed wholly to the iron 
associated with the central enhancement of SN la ejecta. 

Our scheduled Chandra observation of Abell 496 should allow us to trace more 
accurately the gradients in abundances and abundance ratios, given Chandra^s higher 
spatial resolution and relative lack of scattering compared to ASCA. Since our suggested 
mechanism for generating gradients in abundances and abundance ratios should not be 
specific to just one cluster, we will be applying similar analyses to ASCA data for other 
clusters, as well. 

This work was partially supported by the NSF and the State of Alabama through 
EPSCoR grant EIIR-9108761. REW also acknowledges partial support from NASA grant 
NAG 5-2574 and a National Research Council Senior Research Associateship at NASA 
GSFC This research made use of the HEASARC ASCA database and NED. 
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A. Appendix 

We calculate the iron mass to luminosity ratio in the vicinity of the cD and for the 
rest of the cluster. For the cD, we will separately consider the iron from SN la and SN 
II and we will treat the iron from SN la in two ways: we will first attribute to the cD 
all the SN la iron in the central region where the abundance is enhanced; however, this 
may be an overestimate, since the cD's central location makes it difficult to separate 
general intracluster gas from its own interstellar medium, even if the latter is substantial. 
Consequently, we will also consider the possibility that only the "excess" iron at the center 
(that which gives rise to the central abundance enhancement and to the central change in 
elemental abundance ratios) was generated by the cD. This "excess" is one third of all the 
iron within 3', or half of the iron attributable to SN la. 

We derived the intracluster gas density distribution in Abell 496 by using Einstein IPC 
data to determine the shape (core radius and asymptotic slope of a /9-model) of its X-ray 
surface brightness distribution and using the ASCA observations described in §2 to provide 
the flux normalization. Given the iron abundances listed in Table 3, we calculated the iron 
mass within the spherical volume contained between — 3' to be Mpe ~ 5.1 x IO^Mq, as 
listed in Table 5 (also listed are calculations for — 2'). Using our best estimate of the 
SN la iron mass fraction described in §4, we calculated the iron mass from SN la to be 

Fei, ^ 3.3 X lO^Mo 



Mfp, ^ 3.3 X lO^Ma, within - 3'. 



We derived the galaxies' optical luminosity distribution in Abell 496 from the galaxy 
morphological and (V^-band) photometric data of Dressier (1980). We assumed that 
early-type galaxies are the source of the intracluster iron (cf. Arnaud et al. 1992) and 
derived their cumulative spherical luminosity distribution from a deprojection of their 
cumulative surface luminosity distribution. We converted the Dressier (1980) V magnitudes 
to B by assuming B — V = 1 for early-type galaxies. For the central cD, we used the 
photometry of Valentijn (1983), which assigns a total B magnitude of B^ = 13.42, giving 
a blue luminosity of Lb = 2.6 x lO^^h'^^ Lq. We distributed the luminosity of the cD over 
several radial bins, using an r^/^ law with the same effective radius, to avoid an artificial 
luminosity spike at the center. Table 5 lists the iron mass to optical light ratio for SN la 
and II ejecta within various projected regions, 0-2', 0-3' and 3-12'. The errors shown in 
Table 5 include the iron mass and luminosity errors; the fitting errors from the luminosity 
deprojection procedure are relatively small, so are not included. 

It can be seen from Table 5 that the iron mass to light ratio for SN II ejecta is ~2-10 
times smaller at the center than in the outer regions. If ellipticals had protogalactic winds 
driven by SN II, this shows that such a wind in the cD was not suppressed by it being at the 
cluster center (i.e. by being at the bottom of the gravitational potential well of the cluster 
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and in the midst of the highest intracluster gas density). Given that we see no significant 
gradient in SN II ejecta (see §4), this suggests that the SN II ejecta in vicinity of the cD 
is simply the result of a fairly uniformly mixed contamination in the cluster (in the regions 
observed) . 

The nominal SN la iron mass to light ratio in the central region is somewhat less than 
that of the outer parts, but the associated errors are large enough that this difference is not 
significant. Thus, the bulk of the SN la ejecta produced by the cD has been retained in its 
vicinity. 

We also use X-ray observations of an individual elliptical galaxy to see how much 
gaseous iron it has accumulated per unit optical luminosity and compare with the cD in 
Abell 496. NGC 4636 is a particularly X-ray luminous elliptical for its optical luminosity, 
10° from the center of the Virgo cluster in the Virgo Southern Extension (Nolthenius 1993). 
Its metric X-ray to optical luminosity ratio is ~5 times larger than the median for ellipticals 
(White & Davis 1997, 1999), indicating that it has been particularly successful in retaining 
its hot gas. Matsushita et al. (1998) suggest that NGC 4636 may even be at the center 
of a small group of galaxies, with its emission enhanced by group gas. If this is true, it 
makes our conclusion even stronger. The optical and X-ray luminosities of NGC 4636 are 
Lb ~ 3 X 10^^ Lq and Lx ~ 4 x 10^^ erg s^^, adopting a distance of 17 Mpc. Very deep 
ASCA exposures of this galaxy have been analyzed by Matsushita et al. (1997; 1998), who 
found an abundance gradient characterized by a central value of ~0.65 solar (converted to 
the photospheric abundance scale), declining to ~0.2 solar 10' away. 

We used the gas distribution of Matsushita et al. (1998) and the abundance distribution 
of Matsushita et al. (1997) to calculate the iron mass within Tr^s in NGC 4636. This 
encompasses the bulk of its iron content and virtually all of its optical light. We compare 
this to the iron mass and optical luminosity within 3' of the center of Abell 496, which 
encompasses most of the region with enhanced abundances. We find that the SN la iron 
mass to luminosity ratio for the cD is ~40 times greater than in NGC 4636. If we attribute 
to the cD only the "excess" amount of SN la iron at the center compared to the rest of 
the cluster, the SN la iron mass to light ratio in the cD is still ~20 times higher than 
within NGC 4636. The discrepancy is actually even greater than indicated, since we have 
overestimated the SN la iron within NGC 4636 by assuming that all of its iron was produced 
by SN la and we have underestimated the SN la iron in the vicinity of the cD by restricting 
ourselves to within 3', while the SIS data show that abundances are enhanced out to ~ 5'. 
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Figure Captions 

Fig. 1. — a) GIS 2 & 3 spectra for inner (0 — 2') (top) and outer (3 — 12') {bottom) regions; 
solid lines represent the best-fitting isothermal model; individual lines and line complexes 
are identified; b) SIS & 1 spectra for inner (0 — 2') (top) and outer (3 — 12') {bottom) 
regions; solid lines represent the best-fitting isothermal model. 

Fig. 2. — Results of simultaneous isothermal fits to SIS & 1 and GIS 2 & 3 data. Errors 
are 90% confidence limits; 1' ~ 57h^Q. a) Temperature (in keV) distribution; b) Abundance 
(relative to solar) distribution. 

Fig. 3. — Confidence contours for two interesting parameters; abundances and temperatures 
from simultaneous fits of the two innermost (0 — 3' & 3 — 6') regions. The three contours 
correspond to 68%, 90% and 99% confidence limits. 

Fig. 4. — Comparison of derived SN la Fe mass fractions obtained from various observed 
elemental abundance ratios for inner (0 — 2') and outer (3 — 12') regions. The ratios involving 
sulfur have been corrected by reducing the theoretical SN II sulfur production by a factor of 
3.8. The average SN la mass fraction derived in this work is denoted fi. 
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Table 1. Effective Exposure Times 



Spectrometer 


Exposure Time 




(kiloseconds) 


SISO 


29.6 


SIS 1 


23.3 


GIS 2 


39.8 


GIS3 


39.6 



Table 2. Joint Spectral Fits" 



Region ^ 


kT 


Abundance 


x^ 


(arcmin) 


(keV) 


(solar) 




0-3 


q /in+0-06 


0.491^1 


1.07 


0-3^ 


o 7n+0.01 
<5-''-'-0.08 


'^•"-'-'--0.03 


1.06 


3-6 


4 17+0-10 


n S7+0-04 

U-'5'-0.04 


1.04 


6-9 


A /Lfi+0-16 


n qo+0.05 
■^•"^^-0.05 


1.07 


9-12 


A AA+0.26 
4-44_o.24 


r, 49+0.09 


1.02 


6-12 


A /in+0.i3 
4-4U_o.i3 


n q4+0.05 
U-34_o.04 


1.11 


3-12 


A OO + 0.08 

4-^0-0.08 


n •^fi+0-03 


1.07 


0-2 


o 94+0.07 


n ^-^+0-04 

"^•"-''J-0.04 


1.03 


0-2^ 


q q7+0.08 
'^•"^'-0.09 


n f^4+o-05 


1.03 



'^Errors are 90% confidence limits 

Distance from the X-ray center 

^same as above but with an extra coohng 
flow component 



-23- 



Table 3. Individual Elemental Abundances^ 







Region 






Element 


0-2' 


0-3' 


3-12' 


gradient? 





n /IO+0.39 


0.62t°:i 


^7+0-47 




Ne 


n O7+0.45 


o.89i°:i 


7'^+0-48 




Mg 


0.201°:^^ 


0.33t°:^^ 


o.oit°:^f 




Si 


omt'ofr 


0.83l°:ll 


omt'ofs 




S 


n cQ+0.19 


r, 40+0.16 


o.2otnf 


V 


Ar 


<0.3 


<0.3 


o.oitg:^f 




Ca 


r, oq+0.48 

^••-'•J-0.33 


o.22+°:i 


o.3itg:^? 




Fe 


r, c;q+0.05 
^•^•J-0.05 


0.50tHl 


r, qf.+0.03 


V 


Ni 


9 C7+0.67 


2.0lt°JI 


1 04+0-'^'^ 
-"-••J^-o.eo 


V 


^Errors ; 


are 90% confidence limits 
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Table 4. Elemental Abundance Ratios'^ 



Element 
Ratio 



grad 



0-2' 



Region 
0-3' 



12' 



Theory^ 
SN la SN II 



Fe Mass Fraction from SN la 
0-2' 0-3' 3-12' 



0/Fe 




r)qi+0.44 
U.yi_Q44 


1 90+0.39 
^•^-3-0.39 


1 co+0.62 
-^•"-'°-0.60 


Ne/Fe 




1 RA+0-53 


1 77+0.66 
-■-•' '-0.38 


9 nn+0-66 

^••-"^-0.38 


Si/Fe 




1 '^8+0-22 
l.dO_Q22 


1.681^:?^ 


1 09+0-29 
i.y/_Q27 


Ni/Fe 




A oc+1.02 


AMtrrl 


9 00+1.12 


0/Si 




U-t>»_o.29 


n 70+0.58 
U-' -3-0.24 


0-82_o:33 


Si/Ni 


V 


n qq+0.08 
'-'•"-*"-' -0.08 


42+0-1 
'-'•*^-0.09 


fi7+0-27 
U-D'-0.26 


Ne/Si^ 




i.04t°:i? 


1 o^+o-4i 

i-UO_o.24 


1 nc;+0.35 
^-•'-'^-0.31 


S/Fe 


V 


i-iot°i^ 


n Qo+0.22 
0-9^-0.28 


0.54t0;30 


0/S 




n QO+0.42 


1 of;+0.49 
-^•^0-0.53 


2.92ll|« 


Ne/S= 




1 ^0+0.53 


1 01+0.79 
-^•8-L-0.64 


3.72i|?^ 


Si/S^ 


V 


-'-•4o_o.37 


1 79+0.43 
-■■•'^-0.52 


q c;r;+2.03 
•-••"-''-'-1.88 


S/Ni 




n 00+0.06 
"-"•^"^-O.O? 


n 94+0-07 
^•^^-0.08 


o.i9t°:l^ 



0.037 


3.82 


0.006 


2.69 


0.538 


3.53 


4.758 


1.65 


0.068 


1.1 


0.113 


2.14 


0.012 


0.76 


0.585 


2.29 


0.063 


1.67 


0.001 


1.18 


0.919 


1.54 


0.123 


1.39 



0.77^ 



h0.12 
-0.13 

0.39t°iS 
0.651^? 

1 no+o-oo 

1-0.09 



0.87] 



0.86 



n 74+0.09 

"-'•'^-0.07 



0.68^ 



hO.ll 
-0.10 

S4+0-14 

U-<J'l-0.24 
62+0-06 

0.76l°:i 



0.59 
0.26 
0.54 



+0.16 
0.16 
0.14 
0.20 

+0.09 
-0.10 

o.4ol°:l 

0.691^:^^ 
0.4810-19 



'^Errors are propagated la errors 

•^SN la: Nomoto et al (1997a); SN II: Nomoto et al (1997b) 

'^Significantly out of theoretical boundaries 
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Table 5. Gas & Iron Masses from SN la & II 







Region 






0-2' 


0-3' 


3- 12' 


Mpe:. (Mo) 


(1.7 ±0.3) X 109 


(3.3 ± 0.6) X 109 


(1.6 ±0.3) X 10^0 


Mpen (Mq) 


(0.7 ±0.2) X 10^ 


(1.8 ±0.4) X 10^ 


(1.6 ±0.3) X 10^° 


Mpei^/Mpe 


0.70 ± 0.05 


0.64 ±0.05 


0.50 ± 0.06 


ylpei, (solar) 


0.38 ± 0.07 


0.33 ± 0.07 


0.18 ±0.06 


^Fen (solar) 


0.16 ±0.06 


0.18 ±0.06 


0.18 ±0.06 


Mfci./I'Be/so 


0.015 ± 0.009 


0.023 ±0.015 


0.033 ± 0.020 


MFe„/^i?E/so 


0.006 ± 0.004 


0.013 ±0.008 


0.033 ± 0.020 


Mgas (Mq) 


2.2 X 10^2 


5.1 X 10^2 


4.5 X 10^3 


^^E/SO 


1.1 X 10^^ 


1.4 X ll^i 


4.9 X 10^^ 


L Be, so /Mgas 


0.050 


0.028 


0.011 
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